 set scheme s2color
  
 
 

use "Data/Tax Rates/TaxRates", clear

merge 1:1 iso year using "TEMP/Master.dta" , keepusing(direct_real   tax_non_trade_real WAGES)
drop if _merge == 2
drop _merge

encode iso, gen(iso_n)
xtset iso_n year


 gen head_tax_i_real=head_tax_i/WAGES
 gen head_tax_real=head_tax/WAGES
 
 
 
* regressions 

eststo clear  
 
 eststo m1a: reghdfe  tax_non_trade_real head_tax_i if currency == "Shilling", noabsorb  cluster(iso_n year)
eststo m1b:	reghdfe  tax_non_trade_real head_tax_i if currency == "Shilling", absorb(year  ) cluster(iso_n year)

 
eststo m2a: reghdfe  tax_non_trade_real head_tax_i if currency == "Shilling", absorb(year iso_n) cluster(iso_n year)
eststo m2b:	reghdfe  tax_non_trade_real head_tax_i if currency == "Franc", absorb(year iso_n) cluster(iso_n year)
 
 
eststo m3a: reghdfe  tax_non_trade_real  head_tax_i_real ,  absorb(year iso_n) cluster(iso_n year)
eststo m3b: reghdfe  tax_non_trade_real  head_tax_real ,  absorb(year iso_n)  cluster(iso_n year)



 esttab   m1a m1b m2a m2b m3a m3b    	using "Output/Tables/Appendix/APP_TAB_2_2.rtf", replace style(tex) cells(b(star fmt(2)) se(par fmt(2)) P_wald(  fmt(2))) stats(          N r2_a     , fmt(2 2 2   ) labels(  "Observations" "Adjusted  R^2 ")) legend label   starlevels(* 0.1 ** 0.05 *** 0.01) drop(_cons   )
  
 esttab   m1a m1b m2a m2b  m3a m3b    	using "Output/Tables/Appendix/APP_TAB_2_2.tex", replace style(tex) cells(b(star fmt(2)) se(par fmt(2)) P_wald(  fmt(2))) stats(          N r2_a     , fmt(2 2 2   ) labels(  "Observations" "Adjusted  R^2 ")) legend label   starlevels(* 0.1 ** 0.05 *** 0.01) drop(_cons   )



 
